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Recent studies and measurements have shown that the phase stability of the signals 
in the VLF region is very much higher than in the HF spectrum. This fact, along with its 
excellent coverage characteristics, has caused considerable interest in employing this med- 
ium for the wide distribution of standard frequencies and time reference. Basic limitations 
in stability of the received signals are discussed, including path phase distortion, carrier-to- 
noise and envelope delay variations as related to precise synchronization of clocks, and 
highly accurate frequency calibrations. 

Also included is a discussion of the present services of standard frequency and time sig- 
nal stations throughout the world at HF, LF, and VLF. 



1. Introduction 

TIk 1 requirement for better standards for precise 
measurements of frequency and time which are 
readily available has constantly increased with 
advances in various fields of science. The scientific 
unit of time is determined astronomically, and is 
determined to an uncertainty of a few parts in 10 9 
in the course of a year. Since frequency is related 
inversely to time interval, this uncertainty, strictly 
speaking, must be transferred to the specification of 
absolute frequency. Frequency generators, how- 
ever, stable to parts in 10 10 or 10 n now exist, and 
form the basis for relative measurements of frequency 
which for many purposes satisfy the increased re- 
quirements mentioned above. 

The history of the development of our present 
system of time is an extremely fascinating subject 
and it is interesting to note Newton's definition of 
time, "Absolute, true and mathematical time, of 
itself, and by its own nature, flows uniformly on, 
without regard to anything external," [Mach, 1942]. 
Most of man's time systems have been based on an 
attempt to find a reference which "flows uniformly 
on." The very great influence upon our lives and 
surroundings caused by the revolution of the earth 
around the sun and the rotation of the earth on its 
axis has resulted in various definitions 2 of time in 
terms of years (orbital rotation), days (axial rota- 
tion), and subdivision of these units in hours, min- 
utes, and seconds. Time based on apparent solar 
days (from noon to noon) although apparently 
satisfactory as a time reference for many applica- 
tions was shown by astronomical observations to 
have appreciable nonuniformity. 



i Present address of A. D. Watt and R. W. Plush: DECO Electronics, 
Inc., Boulder Division, 8401 Baseline Road, Boulder, Colo. 

2 An interesting and useful description of time and time signals is contained in 
"United States Naval Observatory Circular No. 49," published by the U.S. 
Naval Observatory, Washington 25, B.C. (March 8, 1954). See also "Astro- 
nomical Time," G. M. Clemence, Rev. Mod. Phys. 29, 1 (January 1957) and 
"The Way Things Are," P. W. Bridgman, pp. 135-141 (Harvard University 
Press, Cambridge, Mass., 1959). 



Since apparent solar days arc 4 variable in length 
with season, mean solar time was devised and the 
tabulation of this difference between the apparent 
value and the mean value is called the equation 
of time. This difference has a maximum value a 
little in excess of lb min. Apparent solar days are 
variable 4 in length with season for two reasons: 
partly because of the variation in the angular 
velocity of the sun along' the ecliptic in accordance 
with Kepler's second law of planetary motion, and 
partly because the inclination of the ecliptic to the 
celestial equator introduces a variation in the rate 
at which the projected coordinate of the apparent 
sun moves along the celestial equator. The mean 
sun is a fictitious boctv which moves along the 
celestial equator at a uniform rate equal to the 
mean rate of the apparent sun. Because of the 
fact that a solar time reference (date and time of 
day, etc.) will depend upon the observer's position, 
Universal Time (also known as Greenwich mean 
time) was established with the mean solar time 
reference based on the prime meridian at Greenwich, 
England. The problem of relative position for civil 
purposes has been solved by dividing the earth into 
various time zones usuall3 r differing by one hour 
with lines of demarcation chosen to minimize the 
inconvenience caused by crossing time zones in 
heavily populated areas. 

Three kinds of Universal Times are: UTO, which 
is uncorrected mean solar time; UTl, is Universal 
Time corrected for observed polar motion and repre- 
sents the true angular rotation of the earth about its 
axis; 3 UT2, is Universal Time corrected for both 
observed polar motion and seasonal variation in 
speed of rotation of the earth on an extrapolated 
basis. 

The fundamental unit of time was until very 
recently the mean solar second which was defined as 

„ 4 „r> of a mean solar day. This fundamental unit 
86,400 



3 Because of the increased precision with which the earth's rotation can be 
measured with respect to a reference star rather than the sun, the period of the 
earth's rotation with respect to the vernal equinox, called a sidereal day, is used 
for precise determination of the earth's rotation. 
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of time based upon the mean solar day appeared to 
be entirely satisfactory until astronomical observa- 
tions and advances in the fields of communication, 
electronic frequenc}^ standards, and precision in- 
strumentation showed that in fact the earth was not 
rotating in a strictly constant manner. Even after 
the corrections of UTl and UT2, time based on this 
reference did not flow on uniformly as our definition 
of it would require. A more uniform astronomical 
time based on the yearly motion of the earth about 
the sun, called Ephemeris Time, has long been 
known. The present definition of the second is 

31 556 925 9747 ° f the tro P ical > Tear for Januai T °> 
1900 at 12 hr Ephemeris Time [Markowitz, 1959]. 
This was adopted by the International Committee 
of Weights and Measures in 1956 and ratified by 
the 11th General Conference on Weights and Meas- 
ures in 1960. Discussion of various astronomical 
and atomic times has been given by Markowitz 
[1959]. 

Although quartz crystal clocks were instrumental 
in revealing many of the variations in rotation of 
the earth, it was not until atomic standards were 
placed in use that a new reference of time became 
conceivable. The use of molecular and atomic 
spectral lines as frequency standards [Lyons, 1952] 
has led to the speculation that perhaps here we have 
a reference which can be made highly independent 
of external variable influences and as a result may 
produce a time base of greatly improved uniformity 
and convenience. The difficulty of keeping molec- 
ular and atomic clocks in operation for long periods 
of time has up until recently prevented their use 
as a basic time reference. A determination of the 
atomic transition frequency of cesium in terms of 
Ephemeris Time by Markowitz, Hall, Essen, and 
Parry [1958] indicated a value of 9,192,631,770 
±20 c/s which is the best value available at present. 
The deceleration in the rotation of the earth about 
its axis as determined by the cesium standard is 
in the order of 1.16X10 -5 radians/year 2 , i.e., 0.16 
seconds/year 2 , which is in good agreement with 
values obtained by the moon camera of 0.17 seconds/ 
year 2 for that period (June 1955 to June 1958). 
It is possible [Bullard, 1955] that Ephemeris Time 
and Atomic Time may not have the same rates; 
however, further careful experiments will be re- 
quired to determine this. The stability of the 
best long time interval bases, which at present are 
quartz clocks steered by atomic standards is limited 
to values in the order of one part in 10 10 . Frequency 
comparisons quoted to higher precisions are relative, 
and may be made in two ways. One is to intercom - 
pare two oscillators directly for the short term, 
and the other is to compare one with the mean of 
a group of three or more of the same quality and 
assume that their variations are independent, so 
obtaining a statistical increase in precision for the 
group. The resonance of cesium is known in terms 
of the Ephemeris Second to an accuracy of only 
about 2 parts in 10 9 , the limitation being set by 

4 This means that if an atomic clock and "earth" clock were started together, 
at the end of a year the "earth" clock would be 0.08 sec behind the atomic clock. 



the uncertainty in the determination of the Ephe- 
meris Second. Thus, high quality clocks closely 
steered by atomic standards will not realize the 
presently defined unit of time with any greater 
accuracy than this because of the nature of the 
definition itself. Presently, maximum available 
relative frequency stabilities for periods of several 
hours are in the order of: 1 part in 10 11 for Crystal 
standards and 1 part in 10 12 for the Ammonia Maser. 

2. Basic Limitations of Precision in the Fre- 
quency and Time Standards Via Radio 
Signals 

With the availability of very stable time references 
and frequency standards, it is obvious that con- 
siderable care must be exercised if these standards 
are to be distributed over large areas of the world 
without appreciable deterioration. If the terrestrial 
propagation of radio wave energy was at a constant 
velocity equal to the velocity of light in a vacuum 
and with a noise free background, it would be a 
relatively simple matter to distribute frequency 
and time signals with essentially no loss in stability 
and accuracy. Under these conditions, the received 
frequency would be unaffected and the time signals 
would have a given correctable constant delay of 
T' d =d/v where: d is the distance and v is the 
velocity of light ~3X10 5 km/s. In practice the 
radio energy is transmitted along a path with an 
effective velocity (v) which varies with the charac- 
teristics of and conditions along the transmission 
path so that the actual delay, T d , may vary also. 
It is obvious that the time indicated at a given 
receiving site is T r =T t + T d , where the receiving 
clock can be made to indicate transmitted time T t 
by subtracting T d from the received time T r during 
clock setting. At the high frequencies, 2.5 to 30 
Mc/s, presently employed for the distribution of 
frequency and time signals, multiple propagation 
paths frequently exist with appreciably different 
delays which can introduce large apparent errors 
in the received frequency and time. For any given 
path it may be possible by the proper choice of 
frequencies to have essentially one dominant mode 
of propagation; however, even a single mode may 
have appreciable variations in transmission time. 

2.1. Phase Distortion Limitations 

In figure 1 the straight dashed line represents 
the distortionless transmission in free space where 
the phase delay in radians could be given as 



*'- 



2irfd 

^o 



(1) 



where : / is the frequency in cycles per second, and 
d is the path length expressed in kilometers. 

The actual propagation medium, assuming one 
mode of multiple path to be dominant, may be 
considered as an electrical network with phase and 
amplitude versus frequency characteristics which 
vary as a function of time. The actual time 
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Figure 1. Transmission path phase characteristics. 



delay is 



T a =d/v p 



(2) 



where v v is the average path phase velocity in 
kilometers per second. It is well known in electrical 
circuit theory that various types of network delays 
exist, and if we are to consider the phase delay 
defined in eq (1), v must be expressed in terms of 
path phase velocity rather than envelope velocity. 

The actual path phase velocity relative to the 
velocity of light, v p /v Qj is known to be a function 
of frequency and the path involved. Important 
factors may include ionospheric conditions, ground 
conductivity, and surface roughness. In the VLF 
region, Jean, Taylor, and Wait [1960] have shown 
good agreement between experimentally and theo- 
retically determined values of v p over the 4 to 20 
kc/s frequency range. (See also Wait [1961a] and 
Wait and Spies [1961].) At night they have found 
values of v p /v to be c~1.03 at 4 kc/s, 1.01 at 10 kc/s, 
and 1.003 at 20 kc/s. Many references pertinent 
to the velocity of propagation are contained in 
the paper by Jean, Taylor, and Wait [I960]. In 
general, v p increases with surface conductivity, 
and decreases with increase in ionospheric height 
and earth's surface roughness. 

An actual physical path phase characteristic may 
look somewhat like the solid wavy curve of figure 1 
where a frequency A would be delayed by the phase 
</>!. This phase delay in seconds can be written 
as 

T '=W, (3) 

where T d is the equivalent time delay in seconds. It 
should be pointed out that the particular phase 
versus frequency curve illustrated will only apply 
for a given distance and at a given time and that 
it may vary in its exact position as a function of 
time. The effect of these short term variations 
of phase as a function of time is to limit the precision 
with which a given standard frequency can be 
received. 

To obtain the value of the received frequency 
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Figure 2. Short term propagation time (phase) stability day- 
time paths, normalized to single reflection sky wave. 

from the phase measurements it is clear from the 
definition of frequency that the average departure 
from the transmitted frequency, A/, observed in 
terms of the final and initial phase delays is 
given by: 

(4) 



— ^(T)-«(0) 



2tT 

where T is the time interval, in seconds, over which 
the phase comparison is made. 

The relations governing the transmission p>ath 
phase stability and signal-integration times required 
for specific frequency-comparison precisions have 
been described previously in considerable detail 
[Watt and Plush, 1959]. In general, where possible, 
observations should be made near noon at the 
center of the path. At such time, the phase varia- 
tions, due to diurnal and random short term effects, 
are expected to be minimal. Pierce [1957] has 
experimental evidence over a 5200 km path that 
the total diurnal phase variation of the 16 kc/s 
Rugby signal is surprisingly constant throughout 
the year being in the order of 34 ±1 jusec. Jean 5 
over a 7500 km path has observed that the diurnal 
variation of the phase of the 16 kc/s Rugby signal 
has substantial variations in pattern. The observed 
average phase change of ^42 jusec yields about the 
same effective change in ionospheric height; i.e., 
Ah c^16 km using a method described by Wait 
[1959]. The diurnal patterns observed by Jean 

s A. G. Jean, private communication. See also C. J. Chilton [19611. 
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appear to change systematically with seasons and 
appreciably particularly during ionospherically dis- 
turbed conditions. The daytime transmission delay 
appears to be much less affected by these disturbances 
than the nighttime phase. 

The short time variations in propagation time 
resulting from path phase instabilities can, in 
general, be assumed to be randomly distributed 
about an average time delay T d with a rate of change 
similar to the fade rate. This means that a single 
measured value of the delay ranges roughly between 
T d ± c (T d ), where a (T d ) is the standard deviation 
of the random time delay variable. __ In general, 
there will be different values for both T d and a (T d ) 
for day or night paths. Expected standard devia- 
tions of short time propagation time delays are 
shown in figure 2 as a function of frequency for 
a single reflection daytime sky wave in the VLF and 
LF regions. It should be emphasized that this 
trend does not continue up into the HF region. 
The data are obtained from table B-l [Watt and 
Plush, 1959] with the additional point at 10.2 kc/s 
from Tibbals. 6 

2.2. Time Reference and Clock Setting 

It is obvious that if a single precisely known 
uniform frequency were available at a receiving 
location, that it could be employed in the control 
of a time standard clock, provided: (a) that the 
clock could be properly set initially, and (b) that 
there were no interruptions in the frequency standard 
or clock. Since neither of these can be guaranteed, 
it becomes apparent that some method of establish- 
ing a time or phase reference at the receiving 
location is necessary, and that the problem of trans- 
mitting a precise time reference is much more 

e M. L. Tibbals, private communication . 
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Figure 3. Standard frequency and time signal waveforrns. 



difficult than a standard frequency transmission. 
At the higher frequencies where bandwidth limita- 
tions are not severe, this is accomplished by trans- 
mitting rather short pulses whose leading edges at 
the transmitter are quite steep. The delay of these 
time pulses is also variable but it is not the same as 
the phase delay described in eqs (2) and (3). Since 
the pulses are applied in the form of modulation on 
the rf carrier, the "envelope delay" as defined by 
Nyquist and Brand [1930] in terms of the slope of 
the phase characteristic, d^jdj, is the important 
factor. As we will see in the next section, the HF 
paths are very instable in both phase and envelope 
time delay and the timing accuracy obtainable with 
these pulses is considerably less than that available 
at the transmitter. 

In view of the much more stable phase delay 
conditions in the VLF region as described by 
Pierce [1957] and Allan, Grombie, and Penton [1956] 
along with its excellent coverage- characteristics, 
considerable interest has been exhibited in employing 
this medium for the distribution of standard fre- 
quencies and time reference. The channel band- 
widths available in this frequency region, along with 
the very narrow bandwidths of the transmitting 
antenna, prevent the use of short-time pulses for 
the establishment of coarse time markers which 
could be used in identifying a particular rf cycle. 
An alternate method employing two closely spaced 
and alternately transmitted related frequencies is 
being considered. A description of this method is 
given by Morgan [1961]. The transmitted and 
received waveforms for such a narrow band time 
signal system are shown in figure 3 where it is 
apparent that the envelope of the voltage difference 
produces time reference markers with ambiguities 
spaced by the reciprocal of the frequency difference. 
Coarse time markers may, if desired, be provided 
by the times of frequency switching. On-off keying 
can also be considered for date time code. In a 
dual frequency system, transmitted waveforms are 
given as: 

^i^coso^ (5) 

e 2 =cos oo 2 t (6) 

while the envelope of the phasor sum is given as 



#(0envelope=V2+2 COS (o> 2 — C0 2 )^ 



(7) 

It can be observed that at the transmitter, zero 
time is our phase reference. At the receiver, the 
phases of e\ and e 2 are delayed and these two signals 
are 

e lr =cos [<oi*+0i] (8) 



6 2r = cos [co 2 £+</> 2 ] 



(9) 



Since these two signals can be employed to phase 
lock high quality local oscillators, it is possible to 
have two highly constant amplitude sinusoids as 
indicated. The envelope voltage waveform at the 
receiver is then 

£'(0envelopereceived = V2+2 COS [ (c0 2 - Wl ) t + (<fc-0, ) ] (10) 

where it is evident that the important phase is now 
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4>2—<f>\ and the resulting time delay is 

rp _ 4>2- 4>\ 



2lK/.-/l) 



(ID 



where T b is the envelope delay in seconds. 

If the individual frequencies concerned are sepa- 
rated far enough, the path phase fluctuations 
become independent and the standard deviation of 
T b would be 



y(T b y- 



,2(7(0) 



(12) 



2tt(/ 2 -/ 1 ) 

assuming that cr(</>i) — o-(</> 2 ) = 0- (</>). 

It can be seen that the time variations will decrease 
with an increase in frequency separation. It should 
be emphasized that as the frequency separation is 
reduced the relative phases tend to become corre- 
lated and as a result the envelope reference time 
variation does not increase without limit since the 
phase difference tends to approach zero. Measured 
standard deviations of envelope delay variations 
for a 1 kc/s frequency spacing in the VLF band are 
shown by Casselman, Heritage, and Tibbals [1959] 
to be in the order of ±20 Msec. (See also Stone, 
Markowitz, and Hall, I960.) 7 

Employing the ionospheric roughness parameter 
described by Watt and Plush [1959] the standard 
deviation of short term phase variation; i.e., o- (</>), 
of a VLF carrier received over a 4200 km path is 
found to be in the order of 0.09 radians. When this 
value is used in (12) the expected a(T b ) expressed 
in microseconds is shown in figure 4. Experimental 
values for o-(T b ) obtained from Heritage and 
Tibbals 8 for spacing of 200, 1,000, and 3,000 c/s 
appear to substantiate the general level and trend 
of this curve. 

In order to permit identification of a specific cycle 
of the carrier frequency f h it is obvious that o-(T b ) 
must be less than 0.5/ fi — o-(T d ) which at 20 kc/s will 
be in the order of ±22 /xsec, Until such time as the 
expected reduction in o-{T b ) due to integration for 
time intervals long compared to the fade periods is 
obtained, it appears that frequency separations in 
the order of 1 kc/s will be required for carrier cycle 
identification in the 20 kc/s region. If the fre- 
quency separation chosen is rather large, the 
C/2— /1) markers may be difficult to resolve from the 
on-off keying envelope. This is due in part to the 
normal build up time required by the narrow an- 
tenna bandwidths, see typical shapes in figure 5, 
and the variation in position of this envelope as 
received due to path envelope delay. Should the 
on-off keying phase markers received with time 
delay variations <r(T e ) not be stable enough to re- 
solve C/2— fi) markers, it may be necessary to employ 
an additional frequency providing an intermediate 
reference, viz (/ 3 —j } ) < (J 2 —J x ) . 

Since the precision of frequency comparison is 
defined as the phase jitter divided by the observing 

7 Stone, Markowitz, and Hall in a recent paper [1960] have indicated VLF 
time comparisons of 500 /usee at Washington, D.C., from NBA in Panama. 
Employing a 6 db transmitting bandwidth of 44c/s on fig. 4 and assuming that the 
envelope matching would approach that of two frequencies spaced by the 6 db 
bandwidth we obtain <r(T 6 )^450 ^sec for a path of comparable length. Actually 
the spacing chosen should likely be less which would yield better agreement. 

s Private communication. 



period T(7.3 of Watt and Plush, 1959) we can write 



5=V2*(0)/2Tjfr. 



(13) 



The time delay for a given phase delay is T d =<j)/2irf, 
and since a(T d ) = a(<j>)/2irf } we can write (13) as 



V (T d ) = t T/^2. 



(14) 



Figure 6 shows the manner in which this relationship 
can be employed to determine the variations expected 
in a given clock relative to a standard clock T 
seconds after it has been set to coincidence with the 
standard clock. It is interesting to note that if an 
all daylight radio path such as the 16 kc/s trans- 
atlantic path observed by Pierce [1957] is employed, 
the maximum standard deviation of time difference 
is less than 2 /xsec. 

It should be pointed out that the problem of sub- 
tracting out the mean true transmission time, T d , 
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Figure 6. Clock-time deviations from original setting as a 
function of observing period. ( 8 Pierce 1957). 

whether it be a phase delay or envelope delay, still 
exists. Ionospheric disturbances caused by solar 
related proton storms or meteoric showers have the 
largest apparent effect on the nighttime phase delay 
[Chilton, 1961]. As a result, the midday values of 
T d are expected to be relatively constant from day 
to day with the exception that solar flares effect the 
day values. Over some paths, it is possible that 
there may be appreciable seasonal variations in 
even the midday value of T d . Additional studies 
are required before the errors likely in predicting 
T d can be specified as closely as desired. If two or 
possibly three VLP standard frequency and time 
signal stations are operating simultaneously at 
different rf frequencies but with a common fre- 
quency base as well as a common time reference, it 
is expected that the average delays due to transmis- 
sion time can be obtained with greater accuracy. 
Once the average path delays have been subtracted 
out, a clock at a given receiving site can be set with 
a high degree of precision (probably limited mainly 
by the ability to determine absolute average path 
delay), to the time generated by the network of 
VLF standard frequency and time signal broad- 
casting stations. 

2.3. Carrier to Noise Limitations 

The limitations of noise to precision of frequency 
comparison have been previous^ described [Watt 
and Plush, 1959] and we shall only calculate the 
effects of rms carrier to rms thermal noise ratio on 
time signal accuracy. 

From equation (Cl) of Watt and Plush [1959] 
and (3) 

a(T d )=^- (15) 

(noise) 2y2irfi 



where a(T d ) (noise) is now the standard deviation 
of the time reference in seconds derived from the 
carrier frequency (Jx) phase, and N/C is the rms 
noise to carrier ratio in the receiver. With post 
detection filtering or integration following a linear 
detector, N should be calculated in the effective 
bandwidth of the integrator. If we consider noise 
density or in particular, the noise voltage per 1 kc/s 
bandwidth concept where N\ kc has the dimensions 
volts/ Vkc/s, we obtain 



jB r /l,000 



(noise) 2-yj2wfi C/Ni k c 



(16) 



where B r is the effective receiver bandwidth in 
cycles per second. 

When timing is derived by an envelope beat pat- 
tern between two carriers, the effects of noise are 
given by 



(noise) 



N/G 



M/2-/1) 



or 



„, Vff r /l,000 



(17) 
(18) 



where it is obvious that the time stability for these 
conditions is poorer for a given C/Ni k C /s ratio than 
obtains for a(T d ) (noise). 

We can also determine the stability of a time 
reference obtained from a keyed envelope obtained 
at the one half amplitude point of the leading edge 
with coherent detection. For a single sample this 
is shown in the appendix to be 



1 



(noise, sampled) 



(19) 



where: o-(T e ) is in seconds and B t is the transmitter 
6 db bandwidth in c/s. If the carrier is keyed on 
periodically (for example, once each second) the 
resulting envelope reference can be integrated to 
reduce the effects of noise. The improvement is 
equal to the square root of the number of individual 
pulses, and if the repetition rate is F 



r(T t )- 



1 



{CjNjB^TF 



(noise, integrated) 



or 



<r(T e )=<- 



'Br/lfiOOTF 

(C/N lkc )B t 



(20) 



(21) 



where: jTis the integration period, TF is the number 
of samples averaged, B r is the effective receiver 
bandwidth in c/s. It is apparent from (21) that 
the time reference derived from the envelope 
leading edge improves direct^ as the carrier to 
noise ratio and as the first power of the transmitting 
antenna bandwidth and as the square root of the 
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integration time. In all the derivations involving 
the mis carrier to noise ratios, the significant noise 
background considered is of thermal type. The 
actual background will normally be atmospheric 
noise which in a wide bandwidth has a very different 
statistical nature [Watt and Maxwell, 1957]. The 
effects of atmospheric noise upon radio systems are 
in general quite different from those of thermal 
noise [Watt, Coon, Maxwell, and Plush, 1958]. For 
a pulse sampling system with integration, the re- 
quired mis C/N may be appreciably less than indi- 
cated in eqs. (20) and (21). Hefley 9 [1960] for 
example has found with Loran-C systems that for 
equivalent performance the carrier to atmospheric 
noise ratio can be much less than the carrier to 
thermal noise ratio; however, it must be remembered 
that this system employs a very wide RF band- 
width . 

When carrier phase or frequency stability is con- 
sidered, the narrow bandwidths (0. 1 c/s or less) 
employed will make the atmospheric noise statistics 
the same as thermal noise [Watt and Maxwell, 1957] 
and in these cases the formulas given (15) through 
(18) should apply if linear receivers are employed. 

3. Present Services 

There are at present quite a large number of 
standard-frequency and time-signal stations in use 
throughout the world. VLF, LF, and HF are em- 
ployed for large area coverage with some higher 
frequencies being used for local distribution. A 
description of the stations operating along with 
their location and characteristics is given in Annex 
I of CCIR Report 66 [1956]. 

Annex II of the same document describes the 
characteristics of projected standard frequency and 
time signal stations at four additional locations. In 
general, the high frequency stations all employ pulse 
modulation of the carrier frequency for time signal 
transmissions and many of them in addition have 
standard audiofrequencies employed as modulation 
at various times. The stability as broadcast is 
usually very good, for example with the National 
Bureau of Standards station WWV [1960a] the carrier 
and audiotones have stability of one part in 10 9 at 
all times with normal daily deviations of less than 
two parts in 10 10 . The time signals have essentially 
the same stability and the maximum deviation 
from UT-2 is about ± 50 msec. Final corrections to 
the frequency as broadcast versus the U.S. Frequen- 
cy Standard are published monthly by the National 
Bureau of Standards [1960b] in the Proceedings of 
the IRE while final corrections to the time signals 
versus UT-2 time are published by the U.S. Naval 
Observatory. 

3.1. HF Standard Frequency and Time Services 

Because of the variation in the propagation 
medium at HF, the frequency as received is gener- 
ally much less stable than that transmitted. Ac- 
tually, for high precision measurements it is 
frequently necessary to use long averaging periods 

9 Private communication. 



of 10 to 30 days. Ideally, the following conditions 
should prevail at the time measurements are being 
made; (a) all daylight (or darkness) over radio path, 
(b) no ionospheric disturbance in progress, (c) no 
part of the radio path should pass near either 
auroral zone, (d) single mode propagation should be 
existing. In general, frequency determinations can 
be made either by a direct comparison or by means 
of the time pulses. Some idea of the precision of 
frequency comparison obtainable can be seen from 
figure 7 which shows the precisions available from 
WWV and WWVH over paths ranging from approx- 
imately 2,000 to 8,000 km. Line A shows that for 
short period samples, such as 100 sec, the precision 
available is in the order of 2 parts in 10 7 . If similar 
100-sec samples are made each day, after 30 days 
the precision has improved to one part in 10 8 . It 
should be emphasized that this type of frequency 
comparison does not make full use of the transmis- 
sions and that the improvement with observing 
period is not as great as would be predicted if all 
the actual information received was fully utilized 
for the total periods shown. lane B shows the 
results obtained at Boulder employing time pulses 
compared for approximately 10 sec near noon each 
day and averaged by means of a clock. Here the 
daily averaging time is 24 hr. The point shown for 
30 days is arrived at as follows. Daily comp arisons 
of the local clock versus the received time of WWV 
are used to obtain running 10-day averages of time 
difference. Each day one new daily value is added 
the 1 1-days old value is dropped and the re- 
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sultant 10 time differences averaged. The 30-day 
running averages of the above 10-day averaged 
values are obtained in a similar manner. The pre- 
cisions obtained are appreciably better, being in the 
order of ±2 parts in 10 9 for 1 day and approximately 
±2 parts in 10 l ° for 30 days. 

3.2. LF Standard Frequency Transmissions 

At present there are at least six low frequency 
standard transmissions being operated including 
MSF at 60 kc/s in Rugby, England; WWVB at 
60 kc/s in Boulder, Colo.; DCF77 at 77.5 kc/s in 
Mainningen, Germany; OLP at 50 kc/s in Czecho- 
slovakia; HBJ at 96.04 kc/s in Switzerland; and 
A5XA at 133K kc/s at Fort Monmouth, N.J., 
U.S.A. Other characteristics of these transmissions 
are described in Annex III of CCIR Report 66 
[1956]. Pierce [1957] has shown the standard devia- 
tion of frequency comparison for MSF transmissions 
on a transatlantic path to be in the order of 1 times 
10" 10 for a 30-minute observing period. Further 
calculations for LF transmissions are compared with 
experimental results by Watt and Plush [1959] and 
shown here in figure 8. Pierce [1958] has shown 
average frequency comparison precisions from 
WWVB (60 kc/s) as received at Cruft Laboratory, 
of 1.4 parts in 10 n in an observing period of about 
one hr. From equation (14) this would correspond 
to a short term time variation a (T d ) o^ 0.36 jusec. 

Recent plans to distribute time signals via Loran-C 
navigation system, pulse transmissions at 100 kc/s, 
have been described by Doherty [1961]. Very precise 
time signal distribution in the order of 0.2 /xsec 
is possible within the ground wave coverage range 
of about 1,500 km. By employing sky wave modes ; 
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timing precision in the order of 1 /zsec appears likely 
out to about 3,500 km. 

It is obvious from this and the preceding sections 
that in a limited area the LF systems (50 to 100 
kc/s) will provide for more precise comparison of 
frequency and time than the VLF systems (15 to 
20 kc/s). The propagation loss, however, is much 
higher in the 50 to 100 kc/s band which restricts the 
useful range to distances of about 3,000 to 4,000 
km. LF transmissions at 60 kc/s in England [Steele, 
1955] have shown that timing can be obtained to 
±15 psso, by pulse envelope methods. 

3.3. VLF Standard Frequency Transmissions 

Annex III of CCIR Report 66 [1956] describes 
the transmissions from GBR at 16 kc/s which has 
been used as a basis for much of the pioneer work 
on standard frequency broadcasting in the very low 
frequency region. The precisions available from 
this transmission, and GBZ on 19.6 kc/s, based 
primarily on observations by Pierce [1957, 1958] and 
Allan, Crombie and Penton [1956], are included in 
sections 2.2 and 4. Stabilized transmissions with 
phase locked master-slave relationship from Hawaii 
and San Diego have also served as a valuable source 
of information for VLF path phase stability as 
described by Casselman, Heritage, and Tibbals 
[1959]. Recently interim station broadcasts of the 
U.S. frequency standard have been initiated from 
WWVL on 20 kc/s from near Boulder, Colo. Time 
and stabilized frequency broadcasts from NBA at 
18 kc/s also recently have been initiated. 10 Addi- 
tional VLF transmissions are in the process of being 
stabilized in frequency which should make them 
available for path phase stability measurements. 
In the VLF region both theory and experiment 
show that greater observing times are required for 
a given degree of frequency comparison precision 
than is true in the LF region; however, when long 
range paths are considered the low attenuation rates 
obtaining in the VLF band make emissions in this 
frequency region useful for worldwide coverage. 

4. Characteristics and Expected Coverage 
for a VLF Standard Frequency and Time 
Signal Broadcasting Station 

Based on the previous analysis [Watt and Plush, 
1959] the precision of frequency comparison ex- 
pected as a function of observing period T was 
prepared and is shown as figure 8 [Pierce, Winklei, 
and Corke, I960]. 11 It is obvious when comparing 
figures 7 and 8 that a VLF broadcast will provide 
one or two orders of magnitude improvement over 
present HF broadcasts in precision of frequency 
comparison for a given observing period. This 
great increase in performance has clearly indicated 
the desirability of designing and constructing a 
high quality VLF standard frequency broadcasting- 
station or network of 2 or 3 stations whose primary 

io Precise time and stabilized frequency broadcasts at 18 kc/s over NBA, Sum- 
mit, Canal Zone, were announced in U.S. Naval Observatory Time Service, 
Notice No. 8, 18 November 1959. See also a recent paper by Stone, Markowitz, 
and Hall [I960]. 

n A recent paper by Pierce, Winkler, and Corke [1960] describes results at 16 
kc/s over a 5,200 km path that indicate possible precisions of 2X10 -11 . 
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function is to provide extremely precise standard 
frequency and time signals. 

Because of the desirability of having such a station 
near the U.S. standard of frequency located at 
Boulder, Colo., a study has been made to determine 
design parameters and service obtainable at various 
locations. Seventeen different receiving sites were 
chosen and the total number of measurements obtain- 
able each day for 99 percent of all hours were 
calculated based on a radiated power of 100 kw. 
( Calculations were made for a range of radiated powers 
from 25 to 100 kw and the resulting number of 
measurements obtainable as a function of the cost 
per measurement was determined. The results of 
this study [Watt and Plush, 1959] indicated a 
preference for a radiated power in the order of 100 kw. 

The results of this study are summarized in table 1, 
where the various assumptions made are indicated 
along with a tabulation of the total number of 
measurements obtainable for 99 percent of all hours 
with 100 kw radiated at 20 kc/s 12 for the months 
of June and December assuming a precision of fre- 
quency comparison of 1 part in 10 . The method 
of analysis and some of the assumptions made in 
this analysis are contained in the following sections. 

4.1. Propagation Path Attenuation Rates 
It is relatively well known that the attenuation 
rate (a) to be employed in field strength equations 
such as (2) in the paper by Watt and Plush [1959], 
expressed in decibels per 1,000 km (db/K), is depend- 
ent upon a number of factors including: frequency, 
ionospheric conditions, earth's surface conditions, 
and the earth's magnetic field. 

Experimental and theoretical studies indicate that 
there is a broad minimum in the attenuation- versus- 
frequency curve which is centered somewhere in the 
vicinity of 16 to 18 kc/s. The shape of this attenua- 
tion curve is likely to vary between day and night 
conditions, and is also likely to be different for 
east-to-west and west-to-east propagation. It would 
appear from an analysis of results obtained from 
many different sources that the east-to-west (mag- 
netic) attenuation rate in sonic cases may be as much 
as 1 or 2 db/K greater than the west-to-east attenua- 
tion rates. This effect is expected to be more 
pronounced at night; in addition this difference in 
attenuation rate with direction is expected to be 
greater at the lower frequencies; viz, 10 to 14 kc/s, 
and relatively small in the vicinity of 20 kc/s. These 
effects have been discussed by various authors [H. J. 
Round, T. E. Eckerslev, K. Tremellen, and F. C. 
Lunnon, 1925; K. G. Budden, 1951 and 1952; J. R. 
Wait, 1958 and 1961b; and D. D. Crombie, 1958]. 

Approximate values for the average daytime at- 
tenuation rate (a) in the vicinity of 20 kc/s are: 
sea ^2, average land ^4, estimated, values for perma- 
frost and icecap are 7 and 18. When mixed surface 
conditions are encountered along the path, as is 

i 2 It should be noted that 20 kc/s has been employed in these calculations 
because: (1) 20 kc/s is an internationally assigned standard frequency, (2) the 
power requirements for worldwide coverage are expected to be near a minimum 
at this frequency. (3) considerable propagation data is available in this region. 
Further research on attenuation and phase stability should be conducted over 
the whole VLF spectrum to determine the suitability of other frequencies for 
standard frequency and time distribution. 



usually the case in practice, the total average path 
attenuation will not always be exactly equal to the 
sum of the attenuations expected from each individ- 
ual portion of the path. This may be due in some 
cases to a transfer of energy from the dominant to 
higher order modes caused by the discontinuities in 
conductivity at the surface. It is also likely that 
low conductivity at the earth's surface will produce 
appreciably^ greater attenuation during the day when 
the ionosphere is low than at night when it is higher. 

After determining the various amounts of surface 
types along the individual paths, expected average 
attenuation rates have been assigned for each of the 
receiving locations in table 1. In some cases recipro- 
cal path values have been employed and, in general, 
values quoted are as high as anticipated from a 
consideration of all possible attenuation mechanisms. 

Field strengths anticipated for these particular 
paths have been calculated employing the relation- 
ships given by Watt and Plush [1959]; the results 
are shown in column 5 of table 1 for an assumed 
radiated power of 100 kw at a frequency of 20 kc/s. 

4.2. Expected Availability and Precision of Fre- 
quency Comparison 
The carrier to noise available at the receiving 
location may in some cases be the limiting factor as 
regards the time required to obtain a frequency 
comparison of, say, 1 part in 10 9 . Median atmos- 
pheric noise fields expected for the various receiving 
locations have been obtained from Crichlow, Smith, 
Morton, and Corliss [1955] and the results presented 
in columns 6 and 7 of table 1. Combining the pre- 
vious information regarding carrier field intensities 
with the noise field intensity data, we have obtained 
the minimum expected observation times required 
to compare frequency to 1 part in 10 9 for 100 kilo- 
watts radiated as shown in columns 14 and 15. 

Pierce [1957] has shown that the transmission path 
phase is very stable and that the diurnal variations 
caused by sunrise or sunset along the Rugby to Cam- 
bridge path are surprisingly constant throughout the 
year. During the all sunlight period the path phase 
appears more stable than at night for frequencies of 
16 to 60 kc/s. We have determined, on the basis 
described by Watt and Plush [1959], the minimum 
observing tunes required to obtain a precision of 1 
part in 10 9 based on path phase variations expected 
for both day and night conditions. It should be 
noted that a (0) has been assumed to increase as the 
square root of distance which is what one would 
expect of sections of variable dela}' wave guides 
connected in series. These results are presented in 
columns 12 and 13. 

Although it is obvious that in a case where mini- 
mum observing times based on path phase variabil- 
ity and carrier-to-noise ratio are nearly the same, 
the combined effect will require a longer observation 
time than indicated for either limitation. We have 
neglected this effect and have considered only the 
one condition which requires the longer period. 

In order to get some idea of the relative service 
obtainable from the proposed standard frequency 
and time signal broadcast, the total number of 



625 



o' 

_J 
o 



o 

o 



or 
o 









_i 


tT 


rO 


,J. 


■n}- 


u-> 


O 


rO 


o 


CM 


m 


en 


<* 


Tj- 


en 


m 


O 


on 




00 








CO 


O 


CD 












CO 


CD 


OO 








* 






o 


CM 


^r 
















CD 


rO 


m 




— 


— 






fc 


co Q 
o - 


or 




































a: 
cd 

2- 


as 


oo 


8 


O 

CD 


CM 


K 


CM 


o 


CM 


<X> 


rO 


CO 


CM 


2 




m 


in 


<I 


Si 


en 


O 
CO 


T 


rO 


CM 


H 






% 


v? 


CD 
O 


(X> 


lo 


2 


in 


CM 




S Q 




o 




m 














** 


O 
















2 S 

UJ —1 




2 




— 


















CM 
























rr 


O 


CO 


CO 


en 




CD 


cn 


OJ 


CO 


cn 




rO 






3^ 






ao 


CO 




«J- 


en 


en 


■^r 


CD 




en 






en 
















en 




O 
















o 


ao 














2^ 

2 ^o 




O 




— 


- 














— 


rO 














( 


pj 


m 


en 


^- 


IT) 


-y 


CM 






CM 


o 


m 


m 


^j- 


CO 


rO 


ti- 


CD 




Z 
ID 


31 

cd 


O 


en 












I 


1 


rO 


CM 


rO 




~ 












z 








































en 




o 




tj 


CO 


en 




CD 


en 










O 












ao 


GO 


CO 




CD 


en 




<x> 






en 




















3 






en 














rt 




CC 














o 




en 














































































z 
*■ o 


cc 






































2 ii 


0D 

UJ 


en 

CD 


■* 


_ 


CNJ 


„. 


^ 


CO 


CM 




CD 


rO 


rO 


^ 


cn 








•€> 




O 


O 










GO 


CD 


(0 


o 


o 






rO 


CD 


in 






i ii 


o 


O 


o 


O 


o 


o 


O 


o 


en 


o 


o 


O 


o 


O 


O 




ro 




ti 


o 


O 


o 


o 


o 


o 


o 


o 


o 


CD 


CD 


o 


o 


o 


cD 


CD 


o 




Z CO :* 








































It § 














































































2 QQ 






































a 


x< * — 


z 
5 


<- 

o 


O 


o 


o 


en 
o 


o 


m 


en 
o 


„. 


GO 

en 

8 


CD 

O 
O 


o 


tD 

o 


cp 

o 


cn 
<x> 
O 


t 


r^ 






o 


o 


o 


o 


o 


o 


o 


o 


— 


CD 


O 


CD 


o 


o 


CD 


O 


o 




3E o 






































. 








































f|| 


i 


^_ 




•<r 


co 


^r 










CM 


^ 


ao 




s 


^ 




^ 


rO 


i §s 


o 


o 










rO 


o 




o 










UD 






=E 


o 


o 


o 


o 


o 


d 


o 


CD 


— 


CD 


o 


O 


o 


CD 


O 


CD 


O 




•— «j ^; 








































zm> 














































































i&i 










CO 






























s 


o 


o 


o 


o 












o 


o 


o 












— 


8£g 


o 


o 


o 


o 


o 


o 


o 


o 


in 


o 


o 


o 


CD 


o 


o 





o 






o 


o 


o 


o 


C3 


o 


o 


o 


O 


O 


CD 


o 


O 


CD 


CD 


o 


CD 




s 






































3: 


QC 






































2 


CD 

2 


ao 


en 


QO 


CO 


rO 


o 


— 


rO 


CD 


rO 


ao 


CD 


O 


CO 


-~ 


1^. 


O 




g 

si 


8 


Si 


Si 


2 


= 


OO 


ao 


" 


CD 


CM 


— 


" 


O 


^ 


co 


c5 


* 


en 








































z 






































O 


i 


£ 


CVJ 


^ 


CNJ 


3 


in 


* 


_ 


^. 


o 


_ 


oo 


CD 


o 


ao 


<NJ 


(^ 


CD 




^ 


r— ' 


cb 


NT 


^ 


— 


— 


— 


— 


CM 


"*-' 


ao 


CD 


CM 


CD 


cn 


CNJ 


ao 




J5 




















1 


















JC 


ce 






































<i 

z 
o 


CD 

2 


CXJ 


^_ 


CM 


3 


m 


Tf 


_ 


^ 


O 


_ 


co 


CD 


O 


GO 


<NJ 


^. 


CD 




y 


r-' 


cb 


^ 


o 


— 


— ' 


— 


— 


tJJ 


■«J- 


ao' 


CD 


CM 


CD 


en 


e\i 


co 




Z _l 


S 








































































^ 








































Q 










OO 
































z 
zd 




en 


CO 


•* 




o 






CD 


PO 


ao 


CD 


o 








o 


CO 


<I 


~ 








ao 


CO 


^ 


CX> 


cvi 


5; 


2 


o 


^ 


CD 


O 


in 


en 




-C 






































CD O 


or 

CD 

2 






































d p ^ 


o 


O 


o 


o 


O 


o 


p 


p 


CD 


o 


o 


p 


o 


O 


p 


o 


o 




ez r\ 00 






CO 


o 










O 


CO 


liT 




m 


O 




ao 






c*£^ 




rO 


rO 






K~) 


rO 


rO 


rO 


tJ" 


CM 




rO 


rO 


* 


nO 


* 


* 




§§£ 


Q 










































































O :g "» 


UJ 


o 


O 


o 


o 


CD 


O 


o 


o 


o 


O 


O 


o 


o 


O 


O 


O 


O 


CO 


3 




O 


o 






o 




o 




a 


o 


in 


o 


in 


CD 


m 














ti- 








<s- 


rO 


* 


m 


■«r 






rO 






3= — 






































O 


















































































































d! £ g 










































































U_ ** <3 c 
^ ° -J "5 










































































OZ o UJ ». 






































uj — or =v 


en 








ro 




CD 


o 


o 


en 










*r 








DC QC _ — 

or o ^ 




lO 






rO 


rO 




rO 






CO 


in 


»» 






















































































O Q 






































■< 






































or 




































sz_ J 




































, r 


(ESTIMA 

PATH AT 

DAYLIC 

a,db/IOO 


O 


O 


in 


o 


in 


m 


m 


m 


o 


o 


o 


O 


in 


o 


CD 


in 


m 






" 


in 


" 


5f 


* 


* 


in 


in 


in 


" 


ro 








•vT 


* 


£ 


O 


O 


O 


o 


O 


p 


8 


o 


O 


O 


o 


O 


o 


o 


O 


O 


O 














o 




ao 


m 


ao 


in 




CD 


in 


no 




cS 






rO 








<J- 




o 


oo 




ao 








CD 








in 














en 


CO 










en 




■C- 


co 


cn 




a 






































































































g 






o 








































" 


1 

o 
o 


CD 

o 

z" 
o 
cd 

z 

X 


5 


CD 

Li-T 


1 


ioJ 
>-' 
CD 


co" 


z 


or 

CO 
CO 

ID 

o 
o 


5r 

UJ 
CD 


co" 

z 
en 


J 
| 

UJ 
Q 


* 
O 




rM 
B 
Z 

3 


o" 
z 
<I 

o 


5 
cn 

CD" 

or 

3 

to 

UJ 

z 
<t 


1 
co" 

UJ 

rr 

CO 






CO 


z 


=3 




O 




cr 




*■ 




Z 


z 




* 


or 


3T 


UJ 








o 


3T 








o 


o 






<J 


2 


O 


O 




O 








3= 


3 




*: 


or 


2 




2 


z 


L^- 


CO 






cc 




m 


o 


_ 


CM 


rO 


^ 


in 


CD 


^ 


ao 


en 


o 


_i 


CM 


23 


%f 


m 


CD 


M 




Z 





































8 

jT 

i 

z 



c^ 



or ^1 



£18 



measurements with a precision of 1 part in 10 9 
obtainable per 24 hr for 99 percent of all hours have 
been obtained by dividing the number of hours of 
all-day or all-night over the particular path by the 
appropriate minimum observing times and adding, 
as shown in columns 16 and 17 for the months of 
June and December. It can be seen from this tabula- 
tion that the service provided at various locations is 
appreciably different; however, ic is interesting to 
observe that in all but one case an appreciable 
number of measurements can be made throughout 
each day, and that in the worst case, for transmis- 
sions to New Delhi, at least one comparison with a 
precision of 1 part in 10 9 is expected in each day. 
In many of the closer more favorable locations, it 
would appear that precision of 3 parts in 10 10 or better 
can be expected where observing times of at least 
2 hr or more over daylight paths are possible. The 
maximum precision obtainable that can be expected 
and any locations can readily be calculated using the 
procedures outlined earlier [Watt and Plush, 1959]. 

5, Appendix 

5.1. Timing Accuracy Obtainable With an On-Off 
Keyed Carrier in the Presence of Thermal Noise 

It is well known that the slope of the envelope at 
the one-half amplitude point of a keyed carrier 
passed through a filter with a 6 db bandwidth B t is 



dv 
dt 



^EB t 



(Al) 



where E is the maximum or locked keyed carrier 
voltage, and dvjdt is the slope observed at the one- 
half carrier envelope amplitude point. 

When thermal type noise is present, the envelope 
of the keyed carrier can be expected to have a 
standard deviation in amplitude of a {a) for C/N^>1 
with an envelope detector or for all CjN ratios if a 
synchronous detector is employed. If the point at 
which the noise free keyed carrier crosses the E/2 
point is chosen as a reference T , it can be seen that 
the actual crossing point in the presence of noise will 
vary about this point. The standard deviation of 
this crossing time is defined as <r(T e ), and it is easily 
seen that 



*(a)l<r(T e )=dv/dt=EB t 

a(T e ) = a(a)/EB t 



(A2) 
(A3) 



Since E=^2 G where C is the rms carrier and 
o-(a)=V2 N where iVis the rms noise in the receiver 
IF band pass filter, we can write 



'{T.) = 



1 



(C/N)B t 



(A4) 



When the receiver effective bandwidth is B r (c/s) 
and we use the convention of C/Ni kc , i.e., the rms 
carrier to rms noise in a one kc/s effective band, we 
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obtain for a single sample 



<*{T e ) — 777 



lB r /l,000 



(C/N lkc )B t 



(A5) 
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